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2* (57) Abslracl: The invention relates lo a method for estimation by means of a physical model, of the thermal state of a catalytic 
nitrogen oxide trap, installed in the exhaust system of an internal combustion engine, comprising a modelling of the catalytic trap 
using a chain of perfectly stirred reactors RPAi connected in series and, for a current reactor in said chain, a calculation of the 
temperature of the substrate Ts of the catalytic trap and the temperature of the gases Tgs within and at the outlet from said reactor, 

W using at least the temperature and the mass flow of the gases upstream of the reactor. 
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du substrat Ts du picge cataiytiquc ainsi que de la temperature des gax Tgs au scin ct a la sortie de cc rcactcur, a partir au moins dc 
la temperature des gaz et du debit massique des gaz en amont du rcactcur. 
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Precede et dispositif d' estimation de I'etat thermique d'un piege 
catalytique a oxydes d'azote implante dans la ligne d'echappement 
.d'un moteur a combustion interne, 

5 

L'invention concerne T estimation de r6tat thermique d'un 

piege a oxydes d' azote, ou NOx, implante dans la ligne d'echappement 

d*un moteur a combustion interne, et plus particulierement la 

modelisation du champ thermique longitudinal d'un tel piege ou 

10 catalyseur implant6 dans la ligne d'echappement d'un moteur k 

combustion interne fonctionnant en melange pauvre et dont la 

regeneration est effectuee p6riodiquement. 

Les moteurs k combustion interne produisent des gaz 

d'echappement qui contiennent des substances polluantes 

15 r6glement6es, telles que les oxydes d' azote NOx, les hydrocarbures 

imb rules HC, le monoxyde de carbone CO et les particules qu'il est 

n^cessaire de traiter. 

Les seuils admis pour les 6missions de gas polluants des 

v6hicules automobiles sont reguli^rement revus la baisse dans le 

20 cadre de revolution de r^glementation sur la depollution. 

Pour repondre aux emissions des hydrocarbures imbrules, du 

monoxyde de carbone, des oxydes d' azote et des particules, on place 

dans la ligne d*echappement des systemes de traitement continus, tels 
I 

que des catalyseurs d'oxydation pour traiter les hydrocarbures 
25 imbrules et le monoxyde de carbone, ou pe5riodiques, tels que les 
filtres k particules ou les pifeges catalytiques h oxydes d' azote. 

Ces systemes p^riodiques alternent les phases de stockage des 
polluants et des phases de r6gen6ration des pifeges, c'est-a-dire de 
conversion des polluants stock6s en substances non polluantes. 
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Le piege catalytique a oxydes d* azote, ou plus simplement 
catalyseur, est un monolithe perc6 de canaux ouverts ^ leurs deux 
extremites afin de permettre le passage des gaz d'6chappeinent. 

Ce support est aussi design6 par le terme « substrat » dans la 
5 suite du texte. Les canaux sont ecu verts par utie structure poreuse 
appeUe « washcoat » qui presente une grande surface de contact avec 
les gaz. Sur ce washcoat, sont disposes des in6taux precieux (ex. : 
platine, rhodium) et des metaux alcalins- (ex. : cesium) ou alcalino- 
terreux (ex. : baryum) qui presentent des propri6t6s catalytiques 
10 d' adsorption. 

L' invention vise i resoudre le probleme de T estimation et du 
contrSle de la thermique au sein du catalyseur. En effet, il a 6t6 
observe que ces paramfetres influaient directement sur les efficacit6s 
de stockage et de r6g6n6ration du pifege, sur les quantit6s associ6es 
15 d' oxydes d' azote stock6s ou reduits, et sur le vieillissement thermique 
du substrat, du washcoat et de la phase catalytique. 

L*utilisation de sondes de temperature au sein du pidge h 
oxydes d' azote est relativement compliquee et financi^rement 
rSdhibitoire pour des vehicules produits en grandes series. En effet, au 
20 cott unitaire de chaque sonde se rajoute un cott dfl h T augmentation 
du nombre d' entrees sur le calculateur de contrdle moteur, et un coflt 
pour les Stages mat6riels et logiciels associes. 

De plus, k rheure actuelle, on ne dispose pas de proc6d6 fiable 
(ou peu on6reux) permettant de connaitre les efficacit6s de stockage et 
25 de reduction du piege et les quantites d' oxydes d' azote stockes ou 
r6duits associ6s. La connaissance de la thermique au sein du piege 
catalytique pent permettre de les estimer. 

Par ailleurs, la temperature des gaz en amont du pi&ge n*est pas . 
representative de T^tat thermique du substrat, car en particulier 
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Tinertie thermique du piege et rexothermicite des relations 
thermiques ont une influence non negligeable. 

L'invention a 6galement pour but d'utiliser un estimateur par 
modele physique embarque dans le calculateur de la thermique du 
piege a oxydes d* azote afin d'estimer le vieillissement thermique du 
substrat, du washcoat et de la phase catalytique. 

L'invention a egalement pour but de connaitre la ou les 
temperatures les plus pertinentes pour T estimation des efficacites de 
stockage et de reduction des quantites d' oxydes d' azote associ6es, et 
pour la protection du pifege et T estimation de son vieillissement 
thermique, car la forte capacity calorifique du substrat induit de forts 
gradients thermiques au sein du pi&ge. 

Ainsi, le module thermique selon Tinvention permet de 
connaitre le champ thermique longitudinal du support ou substrat de la 
pliase catalytique du piege k oxydes d' azote. C'est un modMe physique 
unidimensionnel pour lequel on discretise la g6om6trie du pifege 
suivant sa longueur en un faible* nombre de reacteurs parfaitement 
agit6s successifs. 

Dans sa forme la plus simple, le module donne pour chaque 
r^acteur la temperature du substrat h partir de la conuaissance de la 
temperature et du d6bit des gaz en amont du piege. 

Par ailleurs, il simule les ^changes de chaleur par convection, 
prenant ainsi en compte son inertie thermique. 

En variante, les ^changes de chaleur par convection avec Tair 
ambiant autour du pi&ge peuvent etre pris en compte h partir de la 
connaissance de la Vitesse du v6hicule et de la temperature de Fair 
exterieur ou de Tair sous le capot moteur. 

Dans une forme particuliferement elaboree, l'invention pent 
figalement prendre en compte I'exothermicite due aux reactions de 
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conversion des gaz sur la phase catalytique k partir de la connaissance 
des d6bits des substances gazeuses en cause. 

Quelle que soit la variante utilis6e, le mod&le selon T invention 
est un modele en cascade ou chaque reacteur utilise en entr6e la 
5 temperature calcul6e dans le reacteur pr6c6dent. 

En d'autres termes, Tinvention telle que' revendiqu6e propose 
tout d'abord un proced6 d' estimation par module physique embarque 
de r6tat thermique d'un pifege catalytique k oxydes d' azote implant6 
dans la ligne d'6chappement d'un moteur k combustion interne. Ce 
10 proc6d6 comporte une modelisation du pidge catalytique par une 
chatne de r6acteurs parfaitement agit6s, connect6s en s6rie. Par 
ailleurs, pour chaque reacteur courant de la chalne, on calcule la 
temp6rature du substrat du pi&ge catalytique ainsi que la temperature 
des gaz de la sortie de ce r6acteur, k partir au moins de la temperature 
IS des gaz et du debit massique des gaz en amont du rdacteur. 

Le nombre de r^acteurs est de preference inf6rieur k 10, par 
exemple sup^rieur ou egal k 3 et inf^rieur ou 6gal k 6. 

Selon un mode de mise en oeuvre de Tinvention, on simule k 
rinterieur de chaque reacteur les ^changes thermiques par convection 
20 en prenant en compte I'inertie thermique du reacteur. 

On pent 6galement prendre en compte les echanges de chaleur 
par convection avec Tair ambiant autour du piege catalytique k partir 
de la connaissance de la vitesse du vehicule incorporant le moteur et 
de la temperature de Tair exterieur ou de I'air sous le capot moteur. 
25 On pent egalement prendre en compte T influence de la 

variation de la pression des gaz au sein de la ligne d^echappement. 

On pent encore prendre en compte I'exothermicite due aux 
reactions de conversion des gaz sur la phase catalytique, a partir de la 
connaissance des concentrations des substances gazeuses en cause. 



2/5/2007, EAST Version: 2.1.0.14 



wo 2004/092555 PCT/FR2004/000908 



10 



Selon un mode de mise en ceuvre de rinvention, pour le 
r6acteur courant, on calcule les enthalpies de reaction, on estime les 
vitesses des reactions chimiqnes, on resout le systeme d'6quations 
differentielles constitue des bilans matiere pour les especes gazeuses 
et des bilans thermiques pour les gaz et pour le substrat, de fagon a 
obtenir la temp6rature des gaz en sortie du reacteur, cette temperature 
de sortie devenant une variable d'entree pour le r6acteur suivant de la 
chaine. 

L'invention a 6galement pour objet un dispositif d'estimation 
par modele physique embarqu6 de I'Stat thermique d'un piSge 
catalytique i oxydes d* azote implant6 dans la ligne d'echappement 
d'un moteur k combustion interne d'un vehicule, comportant une 
chaJne de reacteurs parfaitement agites connect^s en s6rie mod61isant 
le piege catalytique, et, des moyens de traitement embarques sur le 
IS vehicule et aptes pour un reacteur courant de la chaine, k calculer la 
temperature du substrat du piege catalytique ainsi que la temperature 
des gaz k la sortie de ce reacteur, k partir au moins de la temperature 
des gaz et du d6bit massique des gaz en amont du reacteur. 

Selon un mode de realisation de Tinvention, les moyens de 
20 traitement sont aptes a simuler a I'interieur de chaque reacteur les 
^changes thermiques par convection en prenant en compte I'inertie 
thermique du reacteur. 

Les moyens de traitement peuvent etre figalement aptes k 
prendre en compte les ^changes de chaleur par convection avec Tair 
25 ambiant autour du piege catalytique a partir de la connaissance de la 
Vitesse du vehicule incorporant le moteur et de la temperature de Tair 
ext6rieur ou de Tair sous le capot moteur. 
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Les moyens de traitement peuvent etre dgalement aptes k 
prendre en compte T influence de la variation de la pression des gaz au 
sein de la ligne d'echappement. 

Les moyens de traitement peuvent encore etre aptes a prendre 
5 en compte rexothermicit6 due aux r6actions de la conversion des gaz 
sur la phase catalytique, a partir de la connaissance des concentrations 
des substances gazeuses en cause. 

Selon un mode de realisation de I'invention, les moyens de 
traitement sont aptes, pour le r6acteur courant, k calculer les 
10 enthalpies de reaction, a estimer les vitesses des r6actions chimiques, 
a resoudre le systeme d' Equation differentielle constitue des bilans 
matifere pour les especes gazeuses et des bilans thermiques pour les 
gaz et pour le substrat, de fagon k obtenir la temperature des gaz en 
sortie du reacteur, cette temp6rature de sortie devenant une variable 
IS d'entrde pour le r6acteur suivant de la chaine. 

Le dispositif comporte avantageusement un calculateur 
embarqu6 incorporant les moyens de traitement. 

Le dispositif peut 6tre utilise pour estimer refficacitg de 
stockage et de r6gen6ration du pi^ge ou pour estimer la degradation 
20 des performances de reduction des oxydes d' azote et d'oxydation des 
hydrocarbures et du monoxyde de carbone (effet du vieillissement). 

D'autres avantages et caracteristiques de T invention 
apparaitront k I'examen de la description d^taillee de modes de 
realisation et de mise en oeuvre, nullement limitatifs, et des dessins 
25 annexes, sur lesquels : 

-la figure 1 illustre schematiquement un mode de realisation 
d'un dispositif selon Tinvention ; 

-la figure 2 illustre schematiquement un mode de realisation 
d'un catalyseur selon T invention ; 
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-ia figure 3 illustre schematiquement une discretisation spatiale 
du catalyseur.selon Tinvention ; 

-les figures 4 et 5 illustrent schematiquement- une 
implementation d*un modele embarque selon 1' invention ; et 

-la figure 6 est une table des variables et constantes utilis6es 
dans une estimation par modele selon T invention. 

Sur la figure 1, qui repr6sente un schema du moteur et du 
syst&me de post-^traitement associe, la r6f6rence VH d6signe un 
v6hicule incorporant le moteur ainsi que la ligne d'echappement 
LECH. 

Sur cette ligne d'echappement est dispose un pifege catalytique 
k oxydes d' azote, reference ici NOxTRAP. 

Le moteur est associ6 ici k un turbocompresseur TCM et 
rensemble moteur comporte par ailleurs de fagon classique un 
echangeur de suralimentation ECSA, une vanne VI associ^e k un 
6changeur, un capteur CPTl pour la pression de suralimentation, des 
volets d'air VL. Par ailleurs, une premiere sonde SI comportant un 
debitmetre ainsi qu'une sonde de temperature, est disposee en amont 
du compresseur, tandis qu'une deuxi&me sonde S2 ou estimateur de 
temperature est disposee en amont du pi&ge NOxTRAP. 

L* ensemble est pilots par un calculateur 61ectronique MCL, 
recevant notamment une information RM de regime moteur, et 
incorporant les moyens de traitement selon T invention. 

Comme illustre sur la figure 2, le piege k oxydes d' azote 
NOxTRAP est un monolithe perce de canaux carres ouverts a leurs 
deux extr6mites afin de permettre le passage des gaz d'echappement. 

Ce support est aussi designe par le terme « substrat » dans la 
suite du texte. Les canaux sont converts par .une structure poreuse 
appelee « washcoat » qui presente une grande surface de contact avec 
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les gaz. Sur ce washcoat, sont dispos6s des metaux precieux (ex. : 
platine, rhodium) et des metaux alcalins (ex. : cesium) ou alcalino- 
terreux (ex. : baryum) qui presentent des propri6t6s catalytiques 
d' adsorption. 

Au traitement des oxydes d' azote, vient se rajouter une 
fonction de post-traitement par oxydation du monoxyde de carbone CO 
et des hydrocarbures imbrules HC. 

En exces d'air, les oxydes d* azote 6mis sont oxydes sous 
Taction catalytique du platine qui est Tun des metaux pi6cieux 
formant la phase catalytique, puis fix6s sur des sites actifs d'oxyde de 
baryum sous forme de nitrate de baryum. 

La phase catalytique permet 6galement d'oxyder les 
Itydrocarbures imbrfil^s et les monoxydes de carbone. 

En pr6sence de reducteurs (m61ange riche), les oxydes d' azote 
sont Iib6r6s des sites de baryum et r^duits sous Taction catalytique du 
rhodium qui est Tun des autres mdtaux precieux formant la phase 
catalytique. 

Pour obtenir la forte concentration de reducteurs n6cessaire, il 
faut faire fonctionner le moteur en melange riche. 

L'efficacit6 de stockage en oxydes d' azote du piege est 
fonction de plusieurs paramfetres essentiels, comme la quantit6 de NOx 
stockes NS et la capacity de piegeage NSC du catalyseur. 

La capacity de pi6geage du catalyseur, en gramme par litre du 
catalyseur, repr6sente la quantity admissible pour obtenir une 
efficacite de stockage cumulee pred6finie. Cette capacity est fonction 
de la temperature du substrat mais egalement de la quantite d' oxydes 
d* azote arrivant sur le pifege ainsi que du debit des gaz. 
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L'efficacite de stockage du piege est aussi tres d6pendante de 
rempoisonnement des sites actifs piegeant les oxydes d' azote par le 
soufre contenu dans les gaz d'echappement. 

En effet, les oxydes de soufre' entrant en competition avec les 
5 oxydes d' azote, viennent se fixer sur les sites actifs d'oxyde de 
baryum. Ces oxydes de soufre ne sont pas reduits ou desorbes lors des 
phases de regeneration des oxydes d' azote. Pour purger le piege de ces 
oxydes de soufre, il faut effectuer une regeneration specifique en 
pr6sence de reducteur (milieu riche) et forte thermique. 
10 A partir de la connaissance de la quantity d' oxydes d' azote 

emis par le moteur et de TefficacitS de stockage du catalyseur qui est 
identifiee comme une relation NSC=/(T), on d6duit le remplissage du 
pidge en oxydes d' azote, 

Le catalyseur NOxTRAP se caract6rise aussi par sa capacite k 
15 stdcker de I'oxygdne. Une partie des reducteurs envoyds lors de la 
phase de regeneration r6agit avec Toxyg^ne pi^g^. Cette fraction de 
reducteur ne participe pas ^ la purge des oxydes d' azote du pifege, ce 
qui influe directement sur l'efficacite de la regeneration. 

La capacite k stocker de Toxyg^ne est aussi fonction de la 
20 temperature du pifege et de son vieillissement thermique. 

L'efficacite de regeneration en oxydes d' azote du piege est 
fonction de plusieurs parametres essentiels» tels que la capacite k 
stocker de I'oxygene et la temperature du substrat. 

Les reactions entrant en jeu lors des , regenerations du pi&ge 
25 sont fortement exothermiques. Elles augmentent fortement la 
dynamique et les niveaux thermiques du piege et doivent gtre 
contrSiees. 

II est particulierement avantageux de connaitre Texothermicite 
produite par ces reactions, ce qui permet de mieux apprehender 
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Tefficacite de traitement des oxydes d' azote ainsi que le 
vieillissement du pi^ge. 

En ce qui concerne le vieillissement thermique du piege, il se 
traduit par une diminution irreversible de la surface specifique 
d'echange accessible aux gaz. Le vieillissement thermique impacte 
done le materiau adsorbant et la phase active en rfiduction, modifiant 
ainsi sa capacite k stocker de I'oxyg&ne et des NOx ainsi que sa 
capacite a reduire ces NOx lors des phases de purge. 

invention vise ^ r6soudre le probleme de 1' estimation et du 
contrSle de la thermique au sein du catalyseur, ces parametres influant 
directement sur les efficacites de stockage et de r6gen6ration du pidge, 
sur les quantites associees d' azote stockees ou r^duites, et sur le 
vieillissement thermique du substrat et de sa phase catalytique. 

De plus, la forte capacit6 calorifique du substrat induit de forts 
gradients thermiques au sein du pifege, d'oii la n6cessit6 de connaJtre 
la ou les temperatures les plus pertinentes pour T estimation des 
efficacites de stockage et de reduction, des quantites d* oxydes d' azote 
assocides, et pour la protection du pidge et I'estimation de son 
vieillissement thermique. 

L'inyention propose done de modeliser le catalyseur par un 
nombre reduit de r6acteurs parfaitement agit6s mis en s6rie. 

Un r6acteur parfaitement agite RPAi, tel que repr6sent6 sur la 
figure 3, se caract6rise par le fait que toutes les variables mises en jeu 
dans ce r^acteur sont uniformes (temperature, concentration, .,.). 
Ainsi, la valeur prise en sortie d'un reacteur parfaitement agit6 RPAi 
par une variable, est egale a la valeur de cette variable au sein meme 
du r6acteur parfaitement agite RPAi. 

A titre indicatif, pour un piege NOxTRAP, selon Tinvention, 
dont le volume total pent etre environ de 2 a 3 litres, on pourra utiliser 
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trois a six r€acteurs RPAi, ayant par exemple chacun un volume egal 
au volume total divise par le nombre de reacteurs. 

Ainsi, apres avoir d6coup6 le catalyseur en n elements ou 
r6acteur RPAi, on calcule, pour chaque element i les enthalpies de 
5 reaction et on estime les vitesses des reactions chimiques. 

On resout le systeme d' Equation diff6rentielle constitu6 des 
bilans massiques pour les composes gazeux et pour le substrat. 

On propage ensuite de proche en proche dans le sens de 
r.6coulement des gaz, de sorte que le rfisultat obtenu h la sortie de 
10 r616ment i devient une variable d'entr6e de rel6ment i+1. 

Sous forme de schema-bloc, T implementation du module MM 
selon r invention prend la forme illustr^e schfimatiquement sur la 
figure 4. 

^ Par ailleurs, le sch6ma-bloc de la figure 5 represente sous 
15 forme de module, la succession des calculs r^alis6e au sein de chaque 
r6acteur. 

Sur la figure 4 et la figure 5, la variable Tgs d6signe la 
temperature des gaz en sortie d'un r^acteur, tandis que Tg dfisigne la 
temperature de ces gaz k Tentr^e du r6acteur. 
2® Sur ces sch6mas, et dans la suite du texte, la signification des 

variables et parametres utilises est fournie par la table de la figure 6. 

Par ailleurs, sur les sch6mas des figures 4 et 5, la pression des 
gaz au sein du catalyseur n'apparait pas comme une entr6e, Elle 
pourrait neanmoins figurer comme une variable d' entree du mod&le 
25 dans le cas oil Ton souhaiterait prendre en compte la contre-pression 
effective dans la ligne d'echappement (filtre h particules, ...)• 

On va maintenant decrire plus en detail les hypotheses et 
calculs effectu6s par les moyens de calcul incorpores au sein du 
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calculateur MCL et permettant d*estimer T^tat thermique du 
catalyseur. 

Les hypotheses g6n6rales suivantes sont faites : 
. les transferts de chaleur par conduction et par 'rayonnement 
5 sont neglig6s, 

. la composition molaire des gaz d'echappement est la 
suivante : 

[Nj] = 70%, [CO^] = 10%, [H2O] = 10%, somme des autres 
composants = 10%. 

10 Le corollaire de cette composition est que la masse molaire du 

gaz d'6chappement vaut 29 g/mole, 

• la pression est 6gale k la pression dans la ligne (pression 
atmosph6rique ou contre-pression echappement due k un filtre h 
particules). 

IS . les gaz suivent la loi des gaz parfaits. 

Calculs priliminaires et initialisations 

Ces calculs et initialisations sont effectu^s une seule fois par 
les moyens de traitement et avant Tint^gration des equations 
20 difffirentielles. lis peuvent Stre effectu6s k Tinitialisation du 
calculateur. 



a) Parametres lifis au monolithe 
25 Longueur du monolithe (m) 



V 

mono 

^mn»i/> " "5 • longueur dumonoliflie(m) 

fnono "^mono ^ 
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^mono 



^^mny,n " WTTI * 616inent dc loDgufiur du inonolithe (m) 

mono ^^iUments^catalyseur ^ ' 



Porosite du monolithe (fraction de vide) 

c —1 ^niono 
Peord 

Longueur du coii d'un canal du monolithe (m) 



'^eanattx 



Element de volume du monolithe (m^) 



V 



^Slcments^atedyseur 



Element de surface du monolithe (m^) 

^mono ~ ^ ^ccaial ^canata ^mono 

Element de surface d^^change externe (m^) 



15 b) Initialisations des variables 

II faut initialiser les variables de temperatures ^ la mise sous 
tension du calculateur. On peut par d6faut fixer cette valeur ou utiliser 
la valeur de la temperature d'eau du moteur ou de la temp6rature 
20 exterieure ou encore de la temperature des gaz en amont du piege afin 
de minimiser Terreur d' estimation aprfes la mise sous tension. 
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Aprfes ces calculs preliminaires et I'initialisation, on effectue 
des calculs principaux qui sont maintenant detaill6s. 

Pour chaque rfiacteur, on effectue les bilans et calculs. suivants. 
Les calculs constituent des 6tapes interm6diaires et dependent des 
5 variables r6solues avec les Equations diff^rentielles (bilans 
6nerg6tiqu6s). 

Notations ; 

. La grandeur x(i,t) d6signe la valeur de la grandeur x k I'instant dans 
10 le i'"« RPA. 

. Pour i=0, la temperature Tg(0,t) . correspond k la temperature des gaz 
en amont du piege. 

(1) Bilans Inergltiques 

15 

Equatioii pour la temperature du substrat du catalyseur 

Le terme EXOTHERME qui correspond a un terme source du a la prise 
en compte de la chimie de conversion des composes gazeux, sera 
20 explicit6 plus bas. 

Equation pour la temperature des gaz 

dT^ dT^ 
^vP.Cp, -^dV^p^V^Cp, -^dV-^H,,^, (T, -T,)dS^ 

25 
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(2) Calculs des parametres physiques des gaz 
Masse volumique des gaz d'echappement (kg/m^) 



5 



20 



RgpTg(.i,t) 

Capacity calorifique du gaz d'echappement (J/K/kg) 



CPg(i,t) peut gtre d^fini par interpolation lineaire k partir de la 
temperature et de la pression des gaz. 
10 A la pression atmosph^rique, on peut ^crire la relation suivante : 

Cp,(i,t)=1000+1.510-'(T,_„,,,.,(i.t))+2.210-^(T,_„j,,,.(i.t))^ -lO'^CT, 
(i.t))' -210-»°(T,_,.„,„,(i.t))^ 

15 Viscosite dynamique des gaz d'6chappement (kg/m/s) 

M-gCi.t) peut etre dgfini par interpolation lineaire a partir de la 
temp6rature et de la pression des gaz. 



A la pression atmospherique, on peut 6crire la relation 
suivante : 



Mi.t)=10 *[-0.3+0.0805437Tg(i.t) -6.9523810-'(Tg(i.t))* + 3.61imO-«(T,(i.t))'] 

25 Conductibilite thermique des gaz d'echappement (W/m/K) 

^(i.t) peut etre d6fini par interpolation lineaire k partir de la 
temperature et de la pression des gaz. 
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A la pression atmospheriqiie, on peut ecrire la relation 
suivante ; 

^(i,t)=2,053 10*\Tg(i,t))^'^^^^ -1,996 10"^ 

(3) Calculs des parametres physiques du substrat 

Capacite calorifique du substrat (cordi^rite) (J/K/kg) : 

10 Cpcord(i,t) peut 6tre defini par interpolation lineaire S partir de la 
temp6rature et de la pression des gaz. 

A la pression atmospherique, on peut Ecrire la relation 
suivante : • 



S 



IS 



20 



Cpcord(i>t)=353,8+l,952T,(i,t)-1.5l lO'^ T,^(i,t)+4,863610-^ T,^(i,t) 

(4) Calcul du coefficient d'^change convectif gaz- 
solide 

Vitesse lineaire des gaz dans le catalyseur (m/s) 



25 
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Nombre de Reynolds 



5 Nombre de Prandtl 



10 



Coefficient de transfert gaz/solide a la temperature Tg(J/K/s/m*) 



"'cord TZ~^ ' 



2L. 



'canal 



Cette formule est uniquement valable pour un 6coulement 
laminaire, ce qui est toujours le cas au sein d'un pifege a NOx (section 
15 des canaux ties petite). 



(5) Calcul du coefficient d'echange convectif 
solide-conditions exterieures 

20 



On definit le coefficient d'echange convectif entre le 
catalyseur et I'air ambiant autour de celui-ci par une table 
d'interpolation fonction de la Vitesse du v6hicule. 
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Suivant la position du catalyseur NOxTRAP dans la ligne 
d'6chappenient, on pent annuler on rendre constant ce coefficient. On 
n6glige alors I'influence de la vitesse du vehicule sur I'^change 
convectif vers Textgrieur et on pent alors se passer de cette vitesse 
5 comme parametre d' entree du modele. 

II est 6galement possible de prendre en compte dans le modele 
les degagements de chaleur dus aux reactions chimiques et a la 
conversion des constituants gazeux. 

Pour prendre en compte rexothermicit6 gen6r6e par les 
10 reactions de conversion (reduction/oxydation) des substances gazeuses 
§ur la phase catalytique, on peut modeliser macroscopiquement les 
degagements d'energie a partir du calcul des vitesses de creation ou de 
disparition des especes chimiques, des enthalpies de ces reactions et 
de la connaissance, ou de 1' estimation, des concentrations molaires de 
15 ces substances k rentr6e du catalyseur. 

On va maintenant d6crire les grandes lignes d'une solution 
permettant de prendre en compte cette exothermicite. 

Cette solution se base sur le schema r^actionnel suivant du 
comportement du catalyseur NOxTRAP : 

20 
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(1) oxydation des hydrocarhures 

(2) oxydation du CO 

(3) adsorption de I'oxygine 

(4) desorption de Voxygene 

(5) adsorption des NO2 

(6) desorption des NO2 

(7) reduction des NO2 par 
ritydrogine 

(8) rSduction des NO2 par le CO 

(9) reduction des NO^ pew les 
hydrocarhnres 

(10) oxydation de H^, equivalonte d. 
cells du CO 

(11) oxydation du CO par Voxygene 
de surface 

(12) oxydalioa de Vkydrogene par 
Poxygine de surface 

(13) oxydation des JSC par 
Voxygine de surface 

En utilisant ce schema r6actionnel, on a fait les hypotheses 
suivantes : 

. le schema r6actionnel utilise, a 13 cinetiques, est volontairement 
simplifi6, 

. les reactions interviennent suivant la phase de fonctionnement 
consideree du pi^ge : 

(1), (2), (3), (4), (5), (6), (10) en phase d'adsorption ou 
stockage des NOx, 
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CO+}^02-^C02 . 

0*^r5_,^0^y^02 

2m^N02-^NOl 

iV0|^J1^2SiV+iV<72 

NO2 +2J72 -^y2^2 +2flr2^ 
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CH^ + 30* > H^0 + CO^ + 3S(? 
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(3) . (4), (5). (6), (7). (8), (9), (11). (12), (13) en phase de 
r^g^n^ration, ou purge 

(4) , (6) en phases d'adsorption et de regeneration. 
. On considere qu'il y a 100% de NO, en amont ou au sein du 
catalyseur pour simplifier le schema r6actionnel. 

Pour calculer les concentrations ou teneur molaire en entree du 
catalyseur, on convertit, pour chaque substance, les quantit6s 
massiques estim6es en sortie moteur en concentration volumique ou 
teneur molaire. Le debit massique intervient dans la formula de 
conversion pour convertir des quantites massiques. Pour passer de 
teneur molaire en concentration volumique, la pression et la 
temperature des gaz sont 6galement requises. L'estimation des 
quantites gazeuses est une entr6e du module de temp6rature. Ces 
quantites gazeuses peuvent 6tre obtenues par exemple en utilisant des 
IS tables d' interpolation (en regime, d^bit), 

Au sein du catalyseur, on peut r6soudre les Equations 
diff6rentielles constituees du bilan massique pour toutes les espdces 
r6actives au sein de chaque r^acteur parfaitement agit^. Cette methode 
requiert n6anmoins une charge importante de calcul pour le micrp- 
20 contr61eur MCL. 

Une autre solution consiste k etablir dans une phase de 
calibration, une repartition des concentrations ou teneur molaire des 
substances poUuees en entree du pi^ge dans les differents reacteurs. 
Le calcul des vitesses de creation/disparition des esp&ces r6actives 
25 s'effectue de la maniSre suivante : 

Pour chaque reaction 616mentaire (i) : xA + yB->zC, 

La Vitesse de cr6ation/disparition est definie par la relation : 
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Hypotheses : 

. Les reactions sont 61ementaires 

. Les termes [X] representent les concentrations molaires ou pressions 
partielles des especes X (especes gazeuses ou sites actifs) suivant que 
Ton prenne en compte I'influence de la pression sur les vitesses de 
creation/disparition. 

Les enthalpies de reaction s*ecrivent de cette fagon (rindice j 
correspond au num6ro de la r6action) : 

AHj(i.t)=/(T) 

On d^finit ici par EXOTHBRMB, le terme repr6sentant I'apport 
de chaleur dQ aux reactions de r^duction/oxydation des especes 
contenues dans les gaz sur la phase catalytique du monolithe. II 
correspond k la somme des produits des vitesses de 
cr6ation/disparition par les enthalpies de ces reactions : 



EXOTHERME = 




. On peut limiter le calcul de ce terme en ne prenant en compte 
que les reactions pr6sentant une exothermicit6 pr6pond6rante. On ne 
precise pas ici le choix de ces reactions, 

. Le terme EXOTHERME s'annule dans le cas ott Von ne 
souhaite pas prendre en compte les apports de chaleur dus aux 
r6actions chimiques. 
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REVENDICATIONS 

1- Proc^d6 d*estimation par mod&le physique embarqu6 de Tetat 
thermique d'un piege catalytique k oxydes d' azote implants dans la 
ligne d'echappement d'un moteur k combustion interne, comportant 

5 une modelisation du piege catalytique par une chaine de reacteurs 
parfaitement agites (RPAi) connectes en serie, et, pour un reacteur 
courant de la chaine, un calcul de la temperature du substrat (Ts) du 
pi^ge catalytique ainsi que de la temperature des gaz (Tgs) au sein et a 
la sortie de ce r6acteur, a partir au moins de la temperature des gaz et 
10 du debit massique des gaz en amont du reacteur. 

2- Procede selon la revendication 1, caracterise par le fait que 
le nombre de r6acteurs est inf^rieur a 10, par exemple superieur ou 
6gal ^ 3 et inf^rieur ou egal a 6. 

3.-Proced6 selon la revendication 1 ou 2, caract€ris6 par le fait 
15 qu'on simule a Tinterieur de chaque r6acteur les ^changes thermiques 
par convection en prenant en compte Tinertie thermique du r6acteur. 

4"Proc6d6 selon la revendication 3, caract6ris6 par le fait qu'on 
prend egalement en compte les echanges de chaleur par convection 
avec Tair ambiant autour du pi&ge catalytique k partir de la 
20 connaissance de la vitesse du v6hicule incorporant le moteur et de la 
temperature de Tair ext6rieur ou de I'air sous le capot moteur, 

5-Procede selon Tune des revendications pr6c6dentes, 
caract6ris6 par le fait qu'on prend egalement en compte Tinfluence de 
la variation de la pression des gaz au sein de la ligne d*6chappement. 
2^ 6-Proced6 selon Tune des revendications prec6dentes, 

caract6ris6 par le fait qu'on prend en compte Texothermicitg due aux 
reactions de conversion des composes gazeux sur la phase catalytique, 
k partir de la connaissance des concentrations des substances gazeuses 
en cause. 



2/5/2007, EAST Version: 2.1.0.14 



wo 2004/092555 PCT/FR2 00 4/00090 8 



23 

7- Proc6de selon Tune des revendications pr6cedentes, 
caract6ris6 par le fait que pour le r6acteur courant, on calcule le^s 
enthalpies de reaction, on estirae les vitesses des r6actions chimiques, 
on r6sout le syst&me d'fiquations diff6rentielles constitue des bilans 

5 massiques et thermiques pour les gaz et pour le substrat, de fafon ^ 
obtenir la temp6rature des gaz en sortie du r6acteur, cette temperature 
de sortie devenant une variable d* entree pour le reacteur suivant de la 
chaJne. 

8- Dispositif d' estimation par modMe physique embarque de 
10 I'etat thermique d'un piege catalytique a oxydes d'azote implant6 dans 

la ligne d'6chappement d'un moteur k combustion interne d'un 
v6hicule, comportant une chaine de reacteurs parfaitement agites 
(RPAi) connect6s en serie mod61isant le pi^ge catalytique 
(NOxTRAP), et, des moyens de traitement (MCL) embarques sur le 
15 ve^icule et aptes pour un r6acteur courant de la chaine, k calculer la 
temperature du substrat du pi^ge catalytique ainsi que la temperature 
des gaz au sein et k la sortie de ce reacteur, k partir au moins de la 
temperature des gaz et du debit massique des gaz en amont du 
reacteur. 

20 9-Dispositif selon la revendication 8, caract6ris6 par le fait que 

le nombre de reacteurs est inferieur k 10, par exemple superieur ou 
egal a 3 et inf6rieur ou egal ^6. 

10- Dispositif selon la revendication 8 ou 9, caracterise par le 
fait que les moyens de traitement (MCL) sont aptes a simuler k 

25 I'interieur de chaque reacteur les echanges thermiques par convection 
en prenant en compte Tinertie thermique du reacteur. 

11- Dispositif selon la revendication 10, caract6ris6 par le fait 
que les moyens de traitement (MCL) sont 6galement aptes k prendre en 
compte les echanges de chaleur par convection avec I'air ambiant 

30 autour du pifege catalytique k partir de la connaissance de la Vitesse du 
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vehicule incorporant le moteur et de la temperature de Tair ext^rieur 
ou de Tair sous le capot moteur. 

12- Dispositif selon Tune des revendications 8 k 11, caract6ris6 
par le fait que les moyens de traitement sont egalement aptes k prendre 
en compte T influence de la variation de la pression des gaz au sein de 
la ligne d'6chappement. 

13- Dispositif selon Tune des revendications 8 k 12, caracteris6 
par le fait que les moyens de traitement sont 6galement aptes k prendre 
en compte I'exothermicit^ due aux reactions de conversion des 
composes gazeux sur la phase catalytique, k partir de la connaissance 
des concentrations des substances gazeuses en cause. 

14- Dispositif selon Tune des revendications 8 k 13, caract6ris6 
par le fait que les moyens de traitement sont aptes, pour le reacteur 
courant, k calculer les enthalpies de r6action, a estimer les vitesses des 
reactions chimiques, a rfisoudre le systeme d* equations differentielles 
constitue des bilans massiques et thermiques pour les gaz et pour le 
substrat, de fa9on a obtenir la temperature des gaz en sortie du 
r6acteur, cette temperature de sortie devenant une variable d'entrde 
pour le reacteur suivant de la chaine. 

15- Dispositif selon Tune des revendications 8 k 14, caract6rise 
par le fait qu'il comporte un calculateur embarqu6 (MCL) incorporant 
les moyens de traitement. 

16- Utilisation du dispositif selon Tune des revendications 8 k 
15, pour estimer Tefficacit^ de stockage et de r6g6neration du pidge. 

17- Utilisation du dispositif selon I'une quelconque des 
revendications 8 k 15, pour estimer la degradation des performances de 
reduction des oxydes d' azote et d'oxydation des hydrocarbures et du 
monoxyde de carbone. 
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Symbole 


Nom 


Unit6 


Sigiilflcatiou 


Nature 




pq 




Concentratioos . (ou teneurs molaires) des substances 
contenues dans les gaz en entree du catalyseur 


variables 


Qieord 


CPcord 


J/kR/K 


chaleur specifique apparente de la cordi^rite 


•Table 




CPg 


J/kK/K 


chaleur specifique des gaz d'6chappemeat 


Table 




CTg 


W/m/K 


conductibilite thennique des gaz d'echappement dans un 
Element de catalyseur 


Table 




Dcanau3; 




nombre de canaux par unit6 de surface do monolithe 


constante 




Dclta_Hi 


J/mole/ 


eathalpie de la rdaction i " 


Tables 




dLmono 


m 


eMment de longueur du monolithe 


cnnstnTitft 




dSext 




el6mentde surface d'^cbange ext6rieur 


constante 




dSmono 




element de surface de monolithe 






dVmono 




element de volume de monolithe 




Eaj 


Eai 


J/xnole 


Energie d'activation lice i la reaction i 


V V/ 11 o ul JLL I c 




epcord 


m 


dpaisseur des parois do cordierite 






Epsv 




fraction de vide 






Hcord 


l/nl/KJs 


cocf&cient de transfert eaz / solide 






Hext 


J/s/KJnf 


coefficient d'echanee ext^rieur 


Table 




Lcanal 


m 


longueur du cote d*un canal 


constante 




Lmono 


,m 


longueur du monolithe 




Mgaz 




. 


Masse molairc des i?az d'6chappcment 




Ms 


Mug 


kg/m/a 


viscosite dynamique des gaz d'^chappem'^nt dans un element 
de catalyseur 


Table 




Mr 




nombre d'el^ments de discretisation spatiale 
• 




Pr 


Pr 




nombre dePrandtl 


Variable 




Ptot 


•Pa 


press ion totale des eaz d'echappement 




Qs ^ 


On: 


fi/s 






R, 


Re 




nombre de Reynolds 


■variable 








Constante dss Raz parfaits 




Pcflrd 


rho^cord 


kg/nf 


masse volumzque apparente de la cordierite 


constante 


Ps 




kg/n? 


masse volumique des gaz d'echappement dans un element de 
catalyseur 


Table 


Pmoao 


rhojmoxio 


kR/n? 


masse volumiqus apparente du monolithe 


constante 


1 


q 


mole/irf/ 
s 


Vitesse de la reaction chimique j 


variable 




Smono 


cm? 


surface de la section du monolithe 


constante 


TamontNT 


TamontNT 


K 


Temperature des gaz en amont du catalyseur 


variable 




Text 


K 


Temperature exterieixre 


variable 


Tg. , „ 


Tg^celsius 


^ Celsius 


Temperature des Raz 


variable.^ 


Ts 


Ts 


Kelvin 


Temperature du substrat 


variable 


Vf 


Vf 


m/s 


Vitesse lineaire des gaz dans Is catalyseur 


variable 


Vitesse 


Vitesse 


koi/h 


Vitesse duvShicule 


variable 


V 


Vmoao 


1 ' 


volume du monolithe 


constante 
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